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Abstract 

We investigate the critical behavior of a reaction-diffusion system exhibiting 

a continuous absorbing-state phase transition. The reaction-diffusion system 

strictly conserves the total density of particles, represented as a non-diffusive 

conserved field, and allows an infinite number of absorbing configurations. 

Numerical results show that it belongs to a wide universality class that also 

includes stochastic sandpile models. We derive microscopically the field theory 

representing this universality class. 
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The directed percolation (DP) [p]] universality class is recognized as the canonical example 
of the critical behavior in the transition from an active to a single absorbing state. This 
universality class appears to be very robust with respect to microscopic modifications, and 
non-DP behavior emerges only in the presence of additional symmetries, such as in the case 
of symmetric absorbing states 0, long-range interactions |J, or infinitely many absorbing 
states ||. 

Recently a new universality class of absorbing-state phase transitions (APT) JlJ coupled 
to a non-diffusive conserved field has been identified ||. This class characterizes the criti- 
cal behavior of several models showing absorbing transitions with a dynamics that strictly 
conserves the density of particles, that is represented by a conserved static (non-diffusive) 
field. The models are tuned to criticality by varying the particle density in the initial state, 
and exhibit an infinite number of absorbing states. This universality class is particularly 
interesting because it embraces also the large group of stochastic sandpile models || (and 
in particular the Manna model |7|]) which are the prototypical examples that illustrate the 
ideas of self-organized criticality (SOC) ||. These are driven dissipative models in which 
sand (or energy) is injected into the system and dissipated through the boundaries, leading 
eventually to a stationary state. In the limit of infinitesimally slow external driving, the 
systems approach a critical state characterized by an avalanche-like response. Recently it 
has been pointed out that this critical state is equivalent to the absorbing phase transition 
present in the fixed energy case, that is, in automata that consider the same microscopic 
rules defining the sandpile dynamics, but without driving nor dissipation [pHTT|j. 

The numerical evidence for the existence of such a general universality class || is cor- 
roborated by the observation that all the models analyzed share the same structure and 
basic symmetries; namely, a conserved and static non-critical field dynamically coupled to a 
non-conserved order parameter field, identified as the density of the active particles. These 
observations have led to the conjecture that, in absence of additional symmetries, all stochas- 
tic models with an infinite number of absorbing states in which the order parameter evolution 
is coupled to a non-diffusive conserved field define a unique universality class ||. 
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In this Letter, we study the non-diffusive field limit for the two species reaction-diffusion 
model introduced in Ref. |12[ (see also Ref. ]13[). In this limit the model has a phase transi- 
tion with infinitely many absorbing states and it conserves the total number of particles, that 
is associated to a non-diffusive conserved field. We present extensive numerical simulations 
of the model in two and three dimensions, and determine the full set of critical exponents. 
The obtained values are compatible with the new universality class conjectured in Ref. ||. 
This definitely shows the existence of a very broad universality class that includes reaction- 
diffusion processes, stochastic sandpile models, and lattice gases with the same symmetry 
properties. For the present reaction-diffusion model, it is possible to derive microscopically 
a field theory (FT) description. The resulting action and Langevin equations exhibit the ba- 
sic symmetries that characterize this universality class, and represent the first microscopic 
derivation of a FT for sandpile models. Notably, the resulting FT description recovers a 
phenomenological Langevin approach proposed for stochastic sandpiles [PUTOH- The analysis 
provided here is a very promising path for a coherent description of several nonequilibrium 
critical phenomena now rationalized in a single universality class. 

We consider the two-component reaction-diffusion process identified by the following set 
of reaction equations: 

B -> A with rate k x . (1) 
B + A -> 2B with rate k 2 , (2) 

In this system, B particles diffuse with diffusion rate Db = D, and A particles do not 
diffuse, that is, Da = 0. This corresponds to the limit Da — > of the model introduced in 
Ref. ||12|| . From the rate Eqs. (HD and (0), it is clear that the dynamics conserves the total 



density of particles p = Pa + Pb, where pi is the density of component % = A, B. In this 
conserved reaction-diffusion model, the only dynamics is due to B particles, that we can 
identify as active particles. A particles do not diffuse and cannot generate spontaneously B 
particles. More specifically, A particles can only move via the motion of B particles that 
later on transform into A because of Eq. (0). In absence of B particles, pa can be considered 



a static field. This implies that any configuration devoid of B particles is an absorbing state 
in which the system is trapped forever. 

It is easy to see |12j that the reaction-diffusion process defined by Eqs. ([[]) and (0) 
exhibits a phase transition from an active to an absorbing phase for a non-trivial value of 
the total particle density p = p c , which acts as the control parameter. The critical value p c 
depends upon the reaction rates ki,k2- The nature of this phase transition (whether it is 
first or second order) for Da 7^ appears to be determined by the ratio between Db and Da 



E| ; the static field case (Da = 0), on the other hand, has never been explored. It is clear 
that the static field conserved reaction-diffusion (SFCRD) model allows, for any density p, 
an infinite number (in the thermodynamic limit) of absorbing configurations, in which there 
are no B particles. This is the key difference with respect to the case in which Da 7^ 0, as in 
Ref . |12j . In the latter case a configuration devoid of B particles consists of many diffusing 



A particles, jumping from site to site. In the long run, all particles can visit all sites, and 
therefore, in a statistical sense, all configuration with a fixed number of A's are equivalent 



and the absorbing state can be considered statistically unique . 

The present SFCRD model seems to possess all the required symmetry (stochastic dy- 
namics, many absorbing states, static conserved field) for being part of the universality class 
conjectured in Ref. f|. In order to test this possibility, we have performed numerical simu- 
lations of the model in a <i-dimensional hypercubic lattice with N = L d sites. Each site can 
store any number of A and B particles; that is, our model can be represented by bosonic 
variables. Initial conditions are generated by randomly placing Npa particles A and Np^ 
particles B, corresponding to a particle density p = Pa + Pb '■ The results are independent 
of the particular initial ratio Pa / P^b ^ apart from very early time transients. The dynamics 
proceeds in parallel. Each time step, we update the lattice according to the following rules, 
a) Diffusion: on each lattice site, each B particles moves into a randomly chosen nearest 
neighbor site, b) After all sites have been updated for diffusion, we perform the reactions: 
i) On each lattice site, each B particle is turned into an A particle with probability r\. ii) 
At the same time, each A particle becomes a B particle with probability 1 — (1 — r2)" s , 



where % is the total number of B particles in that site. This corresponds to the average 
probability for an A particle of being involved in the reaction (0) with any of the B particles 
present on the same site. The probabilities r\ and r 2 are related to the reaction rates k\ 
and h% defined in Eqs. (|1|) and (0). The order parameter of the system is ps, measuring 
the density of dynamical entities. For small initial densities p, the system will very likely 
fall into an absorbing configurations with only frozen A particles. For large densities, the 
system reaches a stationary active state with ps ^ 0. 

As we vary p, the system exhibits a continuous transition separating an absorbing phase 
from an active phase at a critical point p c . The order parameter is null for p < p c , and 
follows a power law ps ~ (p — p c )^, for p > p c . The system correlation length £ and time r, 
that define the exponential relaxation of space and time correlation functions, both diverge 
as p — > p c [|TJ. In the critical region the system is characterized by a power law behavior, 
namely £ ~ \p — Pc\~ u± and r ~ \p — p c \ ~"L The dynamical critical exponent is defined 
as r ~ ( z , with z = v\\/p_l- These exponents fully determine the critical behavior of the 
stationary state of the model (see Ref. ||T|). 

We have studied the steady-state properties of the model in d = 2 and 3, by performing 
numerical simulations for systems with size ranging up to L = 512 and L = 125, respectively. 
Averages were performed over 10 4 — 10 5 independent initial configurations. The values 
considered for the rates are r\ = 0.1 and r 2 = 0.5 in d = 2, and r\ = 0.4 and r 2 = 0.5 
in d = 3. From the finite-size scaling analysis for absorbing phase transitions Jl| we obtain 
the critical point (p c = 0.3226(1) in d = 2 and p c = 0.95215(15) in d = 3) and the complete 
set of critical exponents. A detailed presentation of these results will be reported elsewhere. 
In Fig. m we show as an example the order parameter behavior with respect to the control 
parameter A = p — p c , from which it is possible to calculate directly the (3 exponent. The 
results obtained in d — 2 and 3 are reported in Tables I and II and compared with the 
Manna sandpile model in the respective dimension. 

In APT it is possible to obtain more information on the critical state by studying the 
evolution (spread) of activity in systems which start close to an absorbing configuration |15 



In each spreading simulation, a small perturbation is added to an absorbing configuration. 
It is then possible to measure the spatially integrated activity N(t), averaged over all runs, 
and the survival probability P(t) of the activity after t time steps. Only at the critical 
point we have power law behavior for these magnitudes. In the case of many absorbing 
states, the choice of the initial absorbing state is not unique |TjJ. There are several methods 



to perform spreading exponents in this case, and we have followed the technique outlined 
in Ref. ||. This procedure amounts to the study of critical spreading with the so-called 
"natural initial conditions" at p = p c [[0|. The probability distribution P 8 (s) of having a 
spreading event involving s sites, as well as the the quantities N(t) and P(t) can thus be 
measured. At criticality, the only characteristic length is the system size L, and we can 
write the scaling forms P s (s) = s- T °h(s/L D ), N(t) = Pf(t/L z ), and P(t) = t~ s g(t/L z ) JEj . 
The scaling functions f(x), g(x) and h(x) are decreasing exponentially for x ^> 1, and we 
have considered that the spreading characteristic time and size are scaling as L z and L D , 
respectively. In this case simulations were performed for systems of size up to L = 1024 
in d = 2 and L = 200 in d = 3, averaging over at least 5 x 10 6 spreading experiments. 
The new scaling exponents t s ,D,5, and r\ are measured using the now standard moment 
analysis technique [[T7],|I3|], and a consistency check is executed by performing a data collapse 



analysis. The resulting exponents are summarized in Tables I and II. As a further consistency 
check of our results, we have checked that our exponents fulfill all scaling and hyper-scaling 
relations in standard APT. Despite the apparent diversity in the dynamical rules, we can 
safely include the SFCRD model and the Manna model in the same universality class. The 
reported numerical values provide striking evidence for a single universality class, and a 
further check of the conjecture in Ref. ||. 

From a theoretical point of view, the SFCRD allows the microscopic construction of a field 
theory description that will represent also the critical behavior of all models belonging to the 
same universality class. The construction of the FT follows the standard steps outlined in 
the work of Doi, Peliti, and also in Lee and Cardy [|2T|. This technique consists in recasting 



the master equation implicit in Eqs. ([]]) and (§) into a "second quantized form" via a set of 

6 



creation and annihilation bosonic operators for particles A and B on each site. At this point 
it is possible to map the solution of the master equation in a path integral on the density 
fields weighted by the exponential of a functional action S In our case, we can quote 
the elegant results of Ref. ||12|| , just considering that we have Da = 0. The action of the FT 
is thus 

S = J dxdt [i)[d t + (r - DV 2 )}tp + 4>[d t (j) - AV 2 0] 

+ U^ljjty — -0) + U 2 ^((f) + 4>) 

+ vxiP'i) 2 + v 2 ipip(ip(f) - $<f>) + v 3 ijjip<j)<j)} , (3) 

where ip and <fi are auxiliary fields, defined such that their average values coincide with the 
average density of B particles and the total density of particles, respectively, ip and are 
response fields, and the coupling constants are related to the reaction rates fcj. Namely, D 
represents the diffusion coefficient of B particles, A is initially also proportional to D, and 
r is the critical parameter that plays the role of the mass in the FT and is related to the 
difference of the total density with respect to the critical density p c . The critical point in the 
mean-field theory is of course at r c = 0. By standard power-counting analysis, one realizes 
that the reduced couplings Ui/D have critical dimension = 4, while the couplings Vi/D 
have on their part df^ = 2. This means that when applying the renormalization group 
(RG) and performing a perturbative expansion around the critical dimension 4, one could in 
principle drop all the couplings Vi |22j| . The critical parameter of this theory is the density of 



active sites ip, while <fi serves just to propagate interactions. We can exploit some symmetry 
considerations of the FT to relate the physics of the system to the corresponding analytical 
description. In fact, by neglecting irrelevant terms in the power counting analysis, the action 
(0) is invariant under the shift transformation 

4> ^ 4>' = 4> + 6, r — ► r' = r — u 2 S, (4) 

where 5 is any constant. This symmetry has a very intuitive meaning: If we increase 
everywhere the density of the system by an amount 5, we must be closer to the critical point 



by an amount proportional to 5. In other words, this symmetry represents the conserved 
nature of the system. It is also interesting to write the set of corresponding Langevin 
equations (up to irrelevant terms) by integrating out the response fields ^, in the action 
S: 

dtip = D\7 2 ip — rip — U\i\) 2 — U2ip4> + Tf$, (5) 
d t( j) = avV + V<f- (6) 

Here, and r\§ are noise term with zero mean and correlations t)rj^(x', if)) = 

2uiip(x,t)5(x — x')5(t — t'), (r)^(x, fjrj^x', if)) = —U2ip(x,t)6(x — x')S(t — if) and 
{r)^(x , t)rj^(x' , t')) = 0. The noise terms have a multiplicative nature p3 |, that is the stan- 



dard form in APT. Note that the Vi couplings contribute to the noises correlations with 
higher order terms. These equations have a very clear physical interpretation. The field is 
conserved and static, i.e., it only diffuses via the activity of B particles, represented by the 
field ip PH . On its turn, the field ip is locally coupled to the field <f>, but is non-conserved. 
Noticeably, this set of equations recovers the Langevin description that has been proposed 
on a phenomenological level for stochastic sandpiles P,|TU||, with the extra infromation of 
the cross-correlation term (77^77^) . Indeed, the sandpile model has the same basic symme- 
tries of the present reaction-diffusion model, once the local density field p is replaced by the 
local sand-grain (energy) density and the order parameter is identified with the density of 
toppling sites field PJlOfl. It is then natural to expect that the very same basic structure is 
reflected in a unique theoretical description. This observation substantiates the existence of 
a common universality class that embraces stochastic sandpiles, conserved lattice gases and 
reaction diffusion systems with many absorbing states. 

The complete RG analysis of the field theory would allow to extract estimates for the 
critical exponents to compare with simulations in d = 2 and 3. Unfortunately, several severe 
technical problems are encountered in this case. In general, as pointed out in Ref. fT2"|| , the 



couplings Vi become relevant and should be taken into account in the RG analysis. The 
importance of the couplings v , can be argued by the change of the energy shift symmetry 
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form, Eq. in the case of the full action, Eq. (|3|). Second, and more important, is 
the presence of the singular bare propagator for the field 0, that cannot be regularized by 
adding a mass term m 2 00, since it will obviously break the symmetry (f|). This singular 
propagator gives rise to divergences in the RG perturbative expansions, and the results 
of Ref. |T2[ cannot be extended u tout-courf to the limit Da — > 0. In particular, some 
Feynman diagrams in the e-expansion presented in Refs. [|l2|,|n| are proportional to I/Da- 
Hence, the limit Da — * in the theory with Da ^ is non-analytic; any infinitesimal 
amount of diffusion in the energy field renormalizes to a finite value, and definitely changes 
the universality class of the model. Work is in progress to provide a suitable regularization 
that will allow an e-expansion calculation of the FT critical exponents. 

This work has been supported by the European Network under Contract No. ERBFM- 
RXCT980183. We thank D. Dhar, R. Dickman, P. Grassberger, H. J. Hilhorst, M.A. Munoz, 
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TABLES 



Steady state exponents 





P 


P/"± 




z 


Un 


SFCRD 
Manna 


0.65(1) 
0.64(1) 


0.78(2) 
0.78(2) 


0.83(3) 
0.82(3) 


1.55(5) 
1.57(4) 


1.29(8) 
1.29(8) 








Spreading exponents 








r s 


D 


z 


V 


6 


SFCRD 
Manna 


1.27(1) 
1.28(1) 


2.75(1) 
2.76(1) 


1.54(2) 
1.55(1) 


0.29(2) 
0.30(3) 


0.50(2) 
0.48(2) 


TABLE I. Critical exponents for spreading and steady state experiments in d = 2. Fig- 
ures in parenthesis indicate the statistical uncertainty in the last digit. Manna exponents from 
Refs. §011-0. 








Steady state exponents 








P 


P/u ± 


v\_ 


z 


i/i 


SFCRD 
Manna 


0.86(2) 
0.84(2) 


1.39(4) 
1.40(2) 


0.62(3) 
0.60(3) 


1.80(5) 
1.80(5) 


1.12(8) 
1.08(8) 








Spreading exponents 








T s 


D 


z 


V 


5 


SFCRD 
Manna 


1.41(1) 
1.43(1) 


3.32(1) 
3.31(1) 


1.74(2) 
1.75(2) 


0.16(1) 
0.16(2) 


0.76(2) 
0.75(2) 



TABLE II. Critical exponents for spreading and steady state experiments in d = 3. Fig- 
ures in parenthesis indicate the statistical uncertainty in the last digit. Manna exponents from 



Refs. @0JT1 0. 
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FIGURES 
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FIG. 1. Order parameter behavior (stationary density of B particles) as a function of A = p—p c 
for the reaction diffusion model in d = 2 and 3. The slope of the straight lines is (3 = 0.65 in d = 2 
and = 0.86 in d = 3. 
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